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Abstract. When a sandpile relaxes under vibration, it is known that its measured 
angle of repose is bistable in a range of values bounded by a material-dependent 
maximal angle of stability; thus, at the same angle of repose, a sandpile can be 
stationary or avalanching, depending on its history. In the nearly jammed slow 
dynamical regime, sandpile collapse to a zero angle of repose can also occur, as a 
rare event. We claim here that fluctuations of dilatancy (or local density) are the 
key ingredient that can explain such varied phenomena. In this work, we model the 
dynamics of the angle of repose and of the density fluctuations, in the presence of 
external noise, by means of coupled stochastic equations. Among other things, we are 
able to describe sandpile collapse in terms of an activated process, where an effective 
temperature (related to the density as well as to the external vibration intensity) 
competes against the configurational barriers created by the density fluctuations. 
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1. Introduction 

The angle of repose jTj in a sandpile has been a recurrent enigma for physicists, both in 
paradigmatic and in realistic terms. Ideas of self-organised criticality c 2j used the angle 
of repose as a paradigm - it was rather appealing to think that sandpiles relaxed to 
a marginally stable state with a unique critical angle, which resulted from a universal 
and extensive dynamics. Interest in the angle of repose only increased when it became 
obvious to most physicists that this picture was too simplistic [3]; in fact the angle of 
repose is multiply-valued within a certain range, and relaxational behaviour about it 
reveals a great deal about intrinsic length scales within the sandpile jl]. 

We first summarise the relevant phenomenology. Sandpiles formed by the deposition 
of grains on surfaces have sides that are typically inclined at a finite angle to the 
horizontal . This is the angle of repose 9 r : in practice what is seen is that it can take 
a range of values before spontaneous flow occurs, i.e., the sandpile becomes unstable to 
further deposition. The limiting value of the angle before avalanching occurs is known 
as the maximal angle of stability 9 m . The difference between these two angles, which 
is characteristic of a given material pQ, is often referred to as the Bagnold angle [Sj. 
Also, sandpiles show strong hysteresis|| : thus, depending on its conditions of formation, 
a sandpile can either be stable or in motion at any angle 9 such that 9 T < 9 < 6 m . 
This bistable behaviour has been studied theoretically and experimentally |El Ej , but its 
dynamical origin has not been clearly explained so far, to the best of our knowledge. 

Another enigma concerns the dynamics around the angle of repose. Early 
experiments in the physics literature [H] indicated that, when subjected to a large 
vibration intensity, the angle of repose of a sandpile would decay to zero logarithmically 
with time: gently vibrated sandpiles, on the other hand, after an initial logarithmic 
decay, remained 'jammed' at a finite angle of repose for experimentally observable times. 
This led, in a theory [4J that followed, to the suggestion that different mechanisms 
were responsible for decay in the two dynamical regimes. In the strongly vibrated 
sandpile, grains relaxed independently of each other, and the angle of repose decayed 
logarithmically to zero; for gently vibrated sandpiles, there was insufficient inertial 
energy given to individual grains, and collective dynamics were responsible for the 
'jamming' observed at finite angles. The mean- field approach adopted in that work was 
adequate to identify the different rates of decay; a treatment of the fluctuations, needed 
for example to explain bistability or jamming at the angle of repose, was, however, 
beyond its scope. 

In this work, our objectives are, firstly, to propose a possible mechanism for observed 
bistable behaviour; secondly, and equally importantly, to construct a coherent picture, 

§ Sandpiles formed by other processes, e.g. by drainage through a narrow pore within a flat box, also 
have angles of repose which depend on their history; their description can also be included, with minor 
modifications, in our formulation. 

|| This is a result of the 'athermal' nature of sandpiles, i.e., a consequence of the fact that grain sizes 
are too large for the ambient temperature to have any effect on their dynamics. Thus configurations 
that would be dissolved away by Brownian motion in liquids or gases, survive in sandpiles. 
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involving both fast and slow dynamical modes, of relaxation at the angle of repose. 
The present work is based on coupled stochastic equations which are similar, but not 
identical, to those proposed earlier in ||]. The exact analytical treatment to these 
equations will provide, in the slow dynamical regime of interest, a description of the 
dual problems of stability and collapse of the angle of repose. 

We now summarise the main results of the present work. Our basic picture is that 
fluctuations of local density (especially when they occur coherently, in the limit of a 
small perturbation) are the collective excitations responsible for stabilising the angle 
of repose, and for giving it its characteristic width, known as the Bagnold angle 0. 
Density fluctuations can occur because of shape effects [0] or friction [UEHI; they are 
the manifestation in our model of dilatancy, first observed by Reynolds [n] in 1885. 

Our theory below essentially relates the Bagnold angle to the Reynolds' dilatancy. 
In the absence of dilatancy, we suggest that the angle of repose in a vibrated sandpile 
would decay swiftly to zero^f. When dilatancy is present, density fluctuations add 
to the value of the 'bare' angle of repose. In our theory, 'out-of-equilibrium' density 
fluctuations (such as those that might be found at the start of a shaking process) generate 
the maximal angle of stability 6 m ; on the other hand, asymptotic values of density 
fluctuations 'equilibrated' by vibrations at long times, give rise to the the 'typical' angle 
of repose 8 r . The Bagnold angle, which is defined [3] as the difference 



is thus found to be the difference between nonequilibrium and equilibrium values of the 
dilatancy for a given material. 

Consider a sandpile that has relaxed to 9 r in the presence of low noise; it is now 
stabilised by an equilibrium value of the density fluctuations. We ask the question: what 
is the probability, under these circumstances, that the sandpile collapses to a zero angle 
of repose? The configurational landscape we are dealing with is that of grains rather 
close to jamming; density fluctuations involve small (typically intra-cage) displacements 
of grains about their 'equilibrium' positions in this very disordered network. Typical 
configurational barrier heights for collapse under such circumstances would be rather 
high, as they would involve a global rearrangement of grains. This is rather reminiscent 
of the situation close to the 'dynamical transition' in an earlier random graph model 
of granular compaction; in both cases, long-range correlations (corresponding to system- 
wide density fluctuations) need to develop for an appropriate collapse to occur" 1 ". Clearly 
such events would be rare; while they would clearly be facilitated by an increase 
in the external noise, the effect of dilatancy merits more discussion. An increase 
of dilatancy means that density fluctuations are greater; the effect of disorder is 
greater and configurational landscapes, rougher. Since sandpile collapse requires global 

% In liquids, for instance, the absence of dilatancy is the reason why liquid surfaces do not spontaneously 
sustain themselves at a non-zero angle to the horizontal. 

+ In the case of compaction, 'collapse' corresponds to the attainment of the dynamical transition via 
a collapse of excess void space; here it corresponds to the collapse of the sandpile to a zero angle of 
repose. 



56b 



9-. 



0, 




Dynamics at the angle of repose 



4 



rearrangement of grains, it will be less probable in rough (strongly dilatant) landscapes. 
Our analytical results confirm these qualitative observations; we find spontaneous 
collapse to be an activated process, with an effective temperature e that depends on 
the ratio of the external noise to the ambient value of the density fluctuations. 

2. The model 

The dynamics of the angle of repose 6{t) and of the density fluctuations <f>(t) are 
described by the following coupled stochastic equations: 



The parameters a, A 2 are phenomeno logical constants, while r/i(t), r] 2 {t) are two 
independent white noises such that 



Similar equations, which involve spatial rather than temporal coordinates, have 
been written to describe the orientational structure of bridges in granular media; 
these have been extensively analysed in concurrent work ^3] and their results are in 
good agreement with independent simulations. Density fluctuations are key to both 
phenomena; the angle of repose is stabilised by dilatancy, while linear bridges grow at 
the expense of local density fluctuations. 

From a technical viewpoint, the linear equation ()2.2|) for <p(t) is known as an 
Ornstein-Uhlenbeck equation [T4*| UK] , whereas equation (j2.1j) for 9{t) is non-linear, 
as it contains a quadratic coupling to the Ornstein-Uhlenbeck variable <fi{t). In spite 
of this nonlinearity, we have been able to find the equilibrium state of the coupled 
equations (|2.1|) . ()2.2|) analytically - one of the very rare instances where this is possible. 
In the Appendix, we show that the stationary Fokker-Planck equation describing this 
equilibrium state can be solved in closed form by means of a quadratic Ansatz (|A.13j) . 
This ansatz only works for the very special kind of non-linearity here, where (|2.1|) 
involves the square of the Ornstein-Uhlenbeck variable <fi. We have come across a 
similar quadratic Ansatz for a problem with a quadratic coupling between two stochastic 
variables 16. 

We now motivate our equations physically. The first terms in (|2.1j) and ()2.2j) 
are suggested on stability grounds: neither the angle of repose nor the dilatancy are 
allowed to be arbitrarily large for a stable system. The second term in (j2.1)l affirms 
that dilatancy underlies the phenomenon of the angle of repose; in the absence of 
noise, density fluctuations constitute this angle*. The noise in (12.1)1 represents external 
vibration, while that in ()2.2|) is a version of the Edwards compactivity , related as it 
is to purely density-driven effects. In earlier work J4|, these were related via fluctuation- 
dissipation relations to effective temperatures for (decoupled) fast and slow dynamics. In 

* We write a term proportional to <p 2 on symmetry grounds - we would expect it to depend on the 
magnitude rather than the sign of density fluctuations. 



6 = -aO + bcf) 2 + Ai77i(t) 
0= -c0 + A 2V2 {t). 



(2.1) 
(2.2) 



(Vi(t)vAt')} = 28 i:j 8(t-t'). 



(2.3) 
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this work, the inclusion of correlations will be seen to lead to an effective temperature e, 
related now to the ratio of these two noises, in the slow dynamical regime. 

When the material is weakly dilatant (c ^> a), so that density fluctuations decay 
quickly to zero (and hence can be neglected), the angle of repose 9(t) itself obeys an 
Ornstein-Uhlenbeck equation. It relaxes exponentially fast to an equilibrium state, 
whose variance 

2 eq = ^ (2-4) 

(Jj 

is just the zero-dilatancy variance of 9. 

The opposite limit where c < a is of much greater interest. Here, density 
fluctuations are long-lived. When, additionally, Ai is small, the angle of repose has 
a slow dynamics reflective of the slowly evolving density fluctuations. These conditions 
are written more precisely as 

7<1, £<1, (2.5) 

in terms of two dimensionless parameters: 

The second equality follows from (|3.6jl . We see from this that e is essentially the ratio 
between the fluctuations of the angle of repose 9, in the respective limits when it is 
decoupled from, and coupled to, the density fluctuations <p. Giving e the interpretation 
of an effective temperature (see below), we define the regime (|2.5jl as a low-temperature 
and strongly dilatant regime, governed by the slow evolution of density fluctuations. 

3. Dynamics in and out of equilibrium 

Suppose that a sandpile is created in regime (|2.5J) with very large initial values for the 
angle 9 and dilatancy 0o- The initial stage of the dynamics is a transient one; here, 
the noises are negligible, so that the decay is entirely given by the deterministic parts 
of fl2H) and (gH) : 

9(t) = (9 -9 m )e~ at + 9 m e- 2ct , (3.1) 
0(t) = o e- ct , (3.2) 

with 

«U « (3.3) 
a 

Thus, density fluctuations (p(t) relax exponentially, while the trajectory 9(t) has two 
separate modes of relaxation: 

• a fast (inertial) decay in 9(t) ~ 9 e~ at , until 9(t) is of the order of 9 m , 

• a slow (collective) decay in 9{t) w 9 m e~ 2ct . 

When (p{t) and 9{t) are small enough [i.e., (pit) ~ <p eq and 9(t) ~ 9 T , cf. (13. 4j) 
and ([3.6)1 ] for the noises to have an appreciable effect, the above analysis is no longer 
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valid. The system then reaches the equilibrium state of the full non-linear stochastic 
process represented by (j2.1|) and (|2.2|) . An exact analytical investigation of this, for all 
values of the parameters 7 and e, is presented in the Appendix. 

In order to get a feeling for the more qualitative features of the equilibrium state, 
we note first that <p(t) is an Ornstein-Uhlenbeck process with equilibrium variance: 

< = ^- (3-4) 

M c 

We see next that to a good approximation, the angle 9 adapts instantaneously to the 
dynamics of <fi(t) in regime (|2.5j) : 

m « bJ ^- (3.5) 

a 

The two above statements together imply that the distribution of the angle 9{t) is 
approximately that of the square of a Gaussian variable. The typically observed angle 
of repose 8 T is the time-averaged value 

e, = { e U = ^ = ^1. (3.6) 

a ac 

Equation ()3.5|) then reads 

e{t)*e x ^L. (3.7) 

Equation (|3.7|) makes the physics behind the multivalued and history-dependent nature 
of the angle of repose (referred to in the Introduction) rather clear. Its instantaneous 
value depends directly on the instantaneous value of the dilatancy; its maximal (stable) 
value 9 m is noise-independent [cf. (|3.3j) ] and depends only on the maximal value of 
dilatancy that a given material can sustain stablyfj. Sandpiles constructed above this 
will first decay quickly to it; they will then decay more slowly to a 'typical' angle of 
repose 9 r . The ratio of these angles is given by 



e 



i;2 



f = w> (3 - 8) 

so that 9 m ^> 9 r for (f> ^> (p cq . Since spontaneous flow always occurs above 9 m , it is 
known as the angle of maximal stability pQ. 

Below this, i.e., for 9 r < 9 < 9 m , we have a region of bistability which depends 
strongly on sandpile history. The above analysis shows that: 

• Sandpiles submitted to low noise are stable in this range of angles, at least for long 
times ~ 1/c. 

• Sandpiles submitted to high noise [such that the effects of dilatancy become negligible 
in (|2.1j) ] continue to decay rapidly in this range of angles, becoming nearly horizontal at 
short times ~ 1/a. 

This provides rather satisfying agreement with earlier work jB] on tilted sandpiles. 
In that work, model sandpiles were submitted to tilts of varying magnitudes before 

(j Here stability is defined for times of order 1/c, assumed to be extremely long. 



Dynamics at the angle of repose 



7 



being restored to their original (horizontal) state. Large tilting resulted in avalanching, 
while small tilting did not, even when the post-tilt state involved the same angle of 
repose. (The reason for this involved the very different effects in each case on granular 
configurations in the uppermost layers of the sandpile [§].) Bistability at the angle of 
repose was thus clearly manifest. 

Our conclusions are that bistability at the angle of repose is a natural consequence 
of applied noise (tilt jH] or vibration) in granular systems; a sandpile can either be at 
rest, or in motion, at the same angle of repose, depending on its history. 



4. The dynamics of sandpile collapse 

We now examine the probability that under the prolonged effect of low noise, the 
sandpile collapses; thus we look for events where the angle 8(t) vanishes. Such an 
event is expected to be very rare in the regime ()2.5|) : in fact it occurs only if the 
noise rjx(t) in ({23} is sufficiently negative for sufficiently long to compensate for the 
strictly positive term b<p 2 . The fully analytical confirmation of this argument is presented 
in the Appendix. It predicts that the equilibrium probability for 9 to be negative: 

IT = Prob(fl < 0), (4.1) 

scales throughout regime (|2.5j) as 

T(l/4) lU ' 4 s 1 " a^A 1 ( 1 

[see l)A.33|) ]. The scaling function .F(C), given explicitly in (JA.35j) . decays monotonically 
from .F(O) = 1 to JF(oo) = (see Figure |TJ). 

We explore further the regime ( ^> 1, where the 'temperature' e is much lower than 
the 'barrier height' r y 2 . Here, the result (jA.36|) implies that the equilibrium probability 
of collapse vanishes exponentially fast: 

n~expl-^l^l I. (4.3) 



(?)'*) 



2 

Thus complete collapse becomes an activated process: collapse events occur at 
Poissonian times, with an exponentially large characteristic time given by an Arrhe- 
nius law: 

r ~ 1/n ~ exp ^ (7) 7 j • ( 4 - 4 ) 

The stretched exponential in (|4.3|) is intriguing, as it involves a fractional power of 
the usual barrier-height-to-temperature ratio 7 2 /e, and therefore raises questions about 
the nature of the barriers involved. We investigate this further below. 

The equilibrium probability for the magnitude \<p(f)\ of the Ornstein-Uhlenbeck 
process to remain smaller than some 0i -C e q during a time interval T, is known to fall 
off exponentially as 

p^(T)~«p(~§<A (4.5) 
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Figure 1. Plot of the function T(C,) entering the scaling law (|4.2j) of the 
probability II in the regime ()2.5|) . 



We now ask the question: what is the zero-temperature (e = 0) equilibrium probability 
Po{6) for the angle of repose to assume a very small value 9 <C Q-P. Inserting 
4>l q /4>i ~ @r/Q [cf- ()3.7|) ] in the estimate ()4.5|) . with T ~ 1/a (the relaxation time of 
the uncoupled 9 dynamics), we find that this probability is exponentially small: 



P o (0) ~ exp -/c- 







(4.6) 



This somewhat heuristic argument is borne out by the rigorous analysis of the Appendix, 
which leads to k — 1, and also predicts the prefactor [see (jA.30|) . 

In the presence of a low noise intensity (e <C 1), the exponential tail ()4.6|) gets 
convoluted with a narrow Gaussian generated by the noise, whose variance is O'L = e6^. 
This gives 

(6-6^ 

Setting 9 = 0, we get 

0i 



P(9) 



n 



exp 



exp 



2e6l 



d0i. 



(4.7) 



(4.8) 



The saddle-point of the above integral, 

e x « ( 7 £) 1/3 ^ r , (4.9) 

leads to the stretched exponential probability distribution (j4.3jh 

The above suggests a very strong analogy with the famous problem of random 
trapping [IB]- Consider a particle performing Brownian motion in one dimension, with 
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diffusion constant D, amidst a concentration c of Poisson-distributed traps; once a trap 
is reached, the particle ceases to exist. The survival probability S(t) of the particle is 
therefore the probability that it has not encountered a trap until time t. Assuming a 
uniform distribution of starting points, the fall off of this probability can be estimated 
by first computing the probability of finding a large region of length L without traps, 
and then weighing this with the probability that a Brownian particle survives within it 
for a long time t: 

r°° ( 7r 2 Dt\ 
S{t) ~ J exp f -cL - -jy J dL. (4.10) 

The first exponential factor exp(— cL) is the probability that a region of length L is free 
of traps, whereas the second exponential factor is the asymptotic survival probability of 
a Brownian particle in such a regionff. The integral is dominated by a saddle-point at 

whence we recover the well-known estimate 

3 ( n 2 2 ^.\ 1 /3 



S(t) ~ exp ( -- (2iT 2 c 2 Dt) j . (4.12) 

We now elucidate as fully as possible the nature of the analogy, in order to 
get further physical insight into the problem of sandpile collapse. Table ^ shows 
the quantitative correspondence between parameters in both situations. Both the 
survival probability S(t) and the collapse probability II obey the stretched exponential 
laws (|4.12jl and ()4.3jl . with an identical exponent 1/3. Both these anomalous dynamical 
laws are the result of a saddle-point approximation, which represents in each case 
an optimisation procedure. In the trapping problem, large regions without traps are 
improbable, whereas the particle would decay too fast in small ones (it would get 
absorbed at a boundary); the best compromise (|4.11|) is found to scale as L ~ it 1 / 3 . 
In the sandpile problem, angles too far below 9 r are hard to find, as dilatancy would 
resist their existence [cf. (14. 6j) ]: on the other hand, angles that are too large would inhibit 
collapse, given their larger configurational barriers in the face of the noise. Once again, 
the best compromise ()4.9|) scales as 9\ ~ e 1 ^ 3 . 

We use this analogy to develop the following picture for sandpile collapse. Imagine 
that the collapse is visualised as the motion of an effective particle (an exciton, say), 
represented by the collective co-ordinate 0(f). Under the influence of a temperature e the 
exciton diffuses across a rough landscape defined by the <fi excitations; we can consider 
this landscape to be a frozen background, since the decay rate of 9 is much faster than 
that of <p. Valleys are separated by <p barriers whose typical height scales as 7; however, 
sandpile collapse actually involves the traversal of a much lower optimal barrier, given 
by (|4.3jl . The process of sandpile collapse can therefore be visualised as the 9 exciton's 

ff This probability is exp(— Dq 2 i), with q — ir/L coming from the Dirichlet boundary conditions at the 
absorbing endpoints. 
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Table 1. Quantitative correspondence between the derivations of the stretched 
exponential laws (|4.12|) and (|4.3|) in the one-dimensional random trapping and 
sandpile collapse problems. 



search for, and escape across, the rare low barrier (|4.3|) . in a frozen landscape of large 
barriers of typical height 7. 

At a given temperature e, sandpile collapse clearly depends on the nature of the 
density fluctuations. We look at two opposite cases of non-Gaussianness. If density 
fluctuations are peaked around zero (i.e., the material is almost non-dilatant), this 
implies a much flatter, more ordered 0-configurational landscape, easier for the exciton 
to traverse. This would lead to a 'liquid-like' scenario of frequent collapse, where a finite 
angle of repose would be hard to sustain under any circumstances. An explicit example 
is provided by the 7 — > limit, where the collapse probability scales as e 1 / 4 [see (14. 2 j) ]. 

In the opposite case of strong dilatancy (where large values of are more frequent 
than in the Gaussian equilibrium distribution), sandpile collapse is even more strongly 
inhibited. If, for example, \<f)(t)\ is constrained to remain larger than some threshold 
0th, the stretched exponential in ()4.3j) reverts (in the £ < 1 regime considered) to an 
Arrhenius law in its usual form: 

(0th/0cq) 4 ^ ^ 

This would arise in the case of a strongly dilatant material, such as wet sand; angles of 
repose for such materials can be far steeper than usual, and still resist collapse. 

5. Discussion 

In the above, we have looked at a very familiar problem, that of the decay and eventual 
collapse of the angle of repose, using an approach that combines new ideas with very 
traditional concepts such as dilatancy [TT]. Our simple theory suggests that sandpiles 
created at arbitrarily large angles will decay quickly to the maximal angle of stability; 
their subsequent behaviour is bistable, with jamming at a typical finite angle of repose as 
one outcome, or a continuing fast decay to zero, as another. All of this occurs because of 
dynamical competition between the fast dynamics of angle decay and the slow dynamics 
of density fluctuations, especially in the low-noise regime. The collapse of a sandpile in 
the jammed regime is shown to be a rare event; we have obtained exact results for the 
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ensuing activated process, which turns out to have interesting analogies with the well- 
known trapping problem. Using these analogies, we are able to summarise the process 
of collapse as follows: weakly dilatant sandpiles collapse easily, while strongly dilatant 
ones bounce back. 
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Appendix. Equilibrium probability distribution for the process (12. ip . (12.21) 

In this Appendix we derive an exact expression for the double Laplace transform of the 
joint equilibrium probability distribution P(6, <p) of the variables 6 and 0, for arbitrary 
values of the dimensionless parameters e and 7, and especially in the regime ()2.5j) . 

Our starting point consists in writing down the stationary (i.e., time-independent) 
Fokker-Planck equation describing the equilibrium state of the two-dimensional Markov 
process (12. lj) . ()2.2j) . This equation reads [15J 

<)Je +^± = 0, (Al) 



where 



00 



BP BP 



A ^ + A + {ad - 602) f + C0 I + (a + c)p = °- (A - 3) 



(A.4) 



We introduce the dimensionless variables 

t>r (peq 

where 

< = ^, Or = ^ (A.5) 

1 c ac 

have been introduced in (|H.4jl and (jH.fij) . and define the double Laplace transform of 
P(6, 0) as 



L(x, y) = (exp(-x6 - y<f>)) = J J exp(-x6 - y<p) P(0, <p) d6 d0. (A.6) 

In terms of this function, the Fokker-Planck equation ()A.3)) becomes 
(B 2 L BL\ 2 / BL\ 

x (w + ^) =£xL+iy [ yL -^)' (A - 7) 
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n = Prob(# < 0) = / ° P{9) d9 = [ — (A.ll) 



together with the normalisation L(0, 0) = 1. The dimensionless parameters 

c ac 2 A? . . 

t = ? e = TF4' (A ' 8) 

have been introduced in ()2.6j) . 

Our main interest will reside in the distribution of 9, encoded in the Laplace 
transform 

f(x) = L(x,0) = (exp(-x9)), (A.9) 

so that 

P0)= i ^^/(x), (A.10) 

and especially in the probability for 9 to be negative: 

d 

2i7T X 

We first notice that 

L(0,y) = e^ 2 (A.12) 

obeys (jA.7|) . in agreement with the plain observation that the stationary distribution 
of is a Gaussian with unit variance. The full problem can be solved by making the 
Ansatz that L(x,y) keeps the same functional form in y for any fixed value of x, i.e., 
looking for a solution to (|A.7|) of the form 

L{x, y) = f{x) exp (g(x) y^J , (A.13) 

with /(0) = g(0) = 1, where g(x) is the 'x-dependent variance' of (p. Equation (|A.7|) 
boils down to two ordinary differential equations for f(x) and g(x): 

x(g' + 2g 2 ) + 2 1 (g-l) = 0, (A.14) 
f'=(ex-g)f, (A.15) 

justifying thus the validity of the above Ansatz. 

The above equations can be solved as follows. Equation (jA.14j) is a Riccati 
equation ^Hj for g(x), which can be linearised by setting 

= I (t. _ l) 
9 2\ij x) ' 

The new unknown function ip(x) obeys the second-order linear equation 



(A.16) 



whose normalised regular solution reads 
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where /27-1 is the modified Bessel function. Equations (jA.16J) and (jA.15j) then 
respectively yield the explicit expressions 



and 



g{x) = J- v (A.19) 



with 



x J 27 _i(4y / 7x) 

f{x) = h{x)^ 2 l\ (A.20) 



3^ \ 7 / T(2 7 ) (4 7 x) nX 



/«>(*) hnr, = E 



^(27)12,-1(477^)/ ' 1 ' ^ 

The product formula (|A.20j) for f(x) expresses that the variable 9 is the convolution 
of two independent random variables: 

(i) a positive random variable #0 whose distribution, encoded in fo(x), only depends 
on 7; 

(ii) a Gaussian variable with variance e. 

The cumulants c n of 9 can be obtained by means of the series expansion 

We thus get 



Cl = 1 ' C2= 27TT + £ ' C3= (7 + l)(2 7 + l)' 

48(57 + 3) 384(77 + 6) 

c 4 — 7 , ; \ , n , lN0/ri T^TT; c 5 



(7 + 1)(2 7 + l) 2 (2 7 + 3)' (7 + l)( 7 + 2)(2 7 + l) 2 (2 7 + 3)' 

= 3840(42 7 3 + II87 2 + 107 7 + 30) 
C6 ~ (7 + 1)2(7 + 2)(2 7 + 1)3(2 7 + 3)(2 7 + 5)' ( ' ] 

and so on. All the cumulants of 9 coincide with those of 9q, except C2, whose term linear 
in e represents the Gaussian variable in (ii) above. 

Let us first consider the limiting situation 7 = 0, which is reached when the rate c 
vanishes, so that the dynamics of the variable <fi is entirely frozen. The above solution 
simplifies greatly. Indeed only the terms corresponding to n = and n = 1 survive in 
the series representation (jA.18|) for ip{x), so that 

ip(x) = 1 + 2x, (A.24) 

and therefore 

f (x) = (1 + 2x)- l '\ g(x) = -1—. (A.25) 

1 + 2x 

The expression for f (x) shows that #0 is nothing but the square of a normalised Gaussian 
variable, the latter being identified with 0. 
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The probability II for 9 to be negative therefore reads 

U=[^L ^ . (A.26) 
J 2m xy/1 + 2x 

In the regime e <C 1 of most interest, this expression simplifies as 

II « / == « V ; ' r . (A.27) 

V 21tt r(l/4) v 7 

In the general situation where e and 7 are both non-zero, the results (lA~T9l) - (TOT| 
are much more involved. The following situations deserve our attention: 

• For e = 0, i.e., in the absence of noise in (|2.1j) . the variable # reduces to 9q. Its 
probability density reads 

P (9) - f dx ^ ( (47X)7 " 1/2 ^ V2 f A 28) 

PoW -i2i^ e U(27)/ 27 -i(4v^)J ' ( } 



The parenthesis falls off exponentially for x — > +00, confirming thus that 9 is positive. 
As # — > 0, the integral in the right-hand side is dominated by large values of x, where 
the Bessel function can be approximated by its asymptotic expression 

hj-i(z) « -|=, (A.29) 

irrespective of 7. A saddle-point approximation yields the estimate 

(2 7 /£) 7 /8 7 3 \ 1/4 / 7 \ 

Po(0) « Lr, 4r ex P -Z) . (A.30) 



r(2 7 + 1)1/2 9 

The probability density of 9 therefore falls off exponentially fast for 9 <C 7. 
• For £ > 0, the probability for 9 to be negative, 

H - r *L e^ 2 ( (47 " )7 ' 1/2 ^ V2 f A 31) 

hmx* l v r(27)/ 27 . 1 (4 v ^)y' ' ( j 



is non-zero. As e —>■ 0, the integral in the right-hand side is again dominated by large 
values of x, and a saddle-point approximation yields the estimate 

n K JgZC© 1/4 ex P afl!V /3 V (A , 2) 



r(2 7 + 1) 1 / 2 V9W I 2 

The probability II therefore falls off as a stretched exponential for e <C 7 2 . 

• In the 7-^0 limit, the results (jA.30|) and ()A.32|) simplify, as the first factors of both 
expressions go to unity. Furthermore, throughout the regime (|2.5|) . i.e., whenever e 
and 7 are simultaneously small, the probability II obeys a scaling law of the form 

nss fo7S) m c = ^' (A ' 33) 

This result interpolates between (|A.27|) and (|A.32|) . In order to derive the expression 
of the scaling function it is again convenient to come back to the series 
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representation (jA.18j) for ip(x). For 7 < 1, keeping z = finite, that expression 

simplifies as 




Inserting this scaling estimate into (jA.20|) and (jA.ll|) . we obtain after some algebra 



This is a monotonically decreasing function of £, starting from the value .F(O) = 1, so 
that ()A.27|) is recovered. Its fall-off at large values of £, 



agrees with ()A.32|) in the 7 — > limit. 
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